Integral estimator of broadband omnidirectionality 
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By using the notion of wavelength- and angle-averaged reflectance, we assess 
in a systematic way the performance of finite omnidirectional reflectors. We 
put forward how this concept can be employed to optimize omnidirectional 
capabilities. Wc also apply it to give an alternate meaningful characterization 
of the bandwidth of these systems. 
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Floquet-Bloch theory warrants that the behavior of pe- 
riodically stratified media is determined by the trace of 
the transfer matrix of the basic period. Indeed, whenever 
the magnitude of this trace is greater than 2, no waves 
propagate in the structure and then a stop band appears. 

In the context of electromagnetic optics, this no- 
tion is at the basis of photonic crystals^ (that is, one- 
dimensional periodic layered structures), which have 
been attracting a lot of attention because of their 
amazing property of acting as omnidirectional reflectors 
(ODRs): they reflect light at any polarization, any inci- 
dence angle, and over a wide range of wavelengths^ii^. 

Although there are a number of approaches for ensur- 
ing a trace greater than 2 in the basic period^' 9 - ' 10 ' 11 , the 
most feasible design involves two materials with refrac- 
tive indices as different as possible. Such a bilayer system 
is usually designed at quarter-wave thickness (at normal 
incidence), which is enough to guarantee ODR 1 - ' 13 ' 14 . 
However, this assumes perfect periodicity and so requires 
the system to be strictly infinite. Of course, this is 
unattainable in practice and one is led to consider stacks 
of N periods, which are often appropriately called finite 
periodic structures^. One can rightly argue that when 
N is high enough (say, 50 or more), there should be no 
noticeable differences with the ideal infinite casei&. But 
there are commercial ODR designs considering only very 
few periods^ and, in such a situation, the optimization 
of the basic period deserves a careful and in-depth study. 

To shed light on this issue it is essential to quantify the 
ODR performance in a manner that permits unambigu- 
ous comparison between different structures. We hold to 
previous suggestion ;] 18 ' 19 , but to take into due consider- 
ation the key role of the bandwidth, we propose here to 
average the reflectance over all the incidence angles and 
all the wavelengths in the spectral range. With this tool 
at hand, we revisit finite ODRs and characterize their 
properties, addressing, as a side product, a proper pic- 
ture of the omnidirectional bandwidth for these systems. 

We begin by briefly recalling some background con- 



cepts. The basic period of the finite ODR consists of 
a double layer made of materials with refractive indices 
(jil, nn) and thicknesses (d,L, dn), respectively. The ma- 
terial L has a low refractive index, while H is of a high 
refractive index. To characterize the optical response we 
employ the transfer matrix, which can be computed as 
M = Mi Mh, where M/, and M# are the transfer matri- 
ces of each layer, whose standard form can be found in 
any textbook^. 

As we have mentioned before, band gaps appear when- 
ever the trace of the basic period satisfies 



Tr(M)| > 2. 



(1) 



This condition should be worked out for both basic po- 
larizations. However, it is known that whenever Eq. ([1]) 
is fulfilled for p polarization, it is always true also for 
s polarization. The p-polarization bands are more strin- 
gent than the corresponding s-polarization ones^ and we 
thus restrict our study to the former. 

We next consider an iV-period structure whose basic 
cell is precisely the bilayer LH. We denote this as [LH] N 
and its overall transfer matrix is simply M N , from which 
calculating its reflectance is straightforward. Ac- 

cording to our previous discussion, we average over all 
the incidence angles and all the wavelengths in the spec- 



tral interval AA = A,, 



Amin of interest 
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and take this as an appropriate figure of merit to assess 
the performance as an ODR. Once the materials have 
been chosen, 1Z is a function exclusively of the layer 
thicknesses. 

As a case study, we take the materials to be cryolite 
(NasAlFg) and zinc selenide (ZnSe), with refractive in- 
dices riL = 1.34 and Uh = 2.568, respectively, at the 
wavelength A = 0.65 /im. The spectral window consid- 



ered is from A n 



0.5 /xm to A n 



0.8 /im. In this 
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Fig. 1. (Color on line) Averaged reflectance 1Z in 
terms of number of periods for the optimized [LH] N sys- 
tems (blue triangles) and for the same system when all 
the thicknesses can be different (red circles). Solid lines 
represent the fittings to these points. 



range, the refractive index of the cryolite can be con- 
sidered, to a good approximation, as constant, while for 
the zinc selenide we use the Sellmeier dispersion equation 
n 2 H (X) = 4 + 1.9A 2 /[A 2 - (0.336) 2 ], where A is expressed 
in microns. 

Many commercial packages are available to perform 
layer optimization. In the case of ODR, common meth- 
ods optimize a merit function that (quadratically) mea- 
sures how the calculated reflectance separates from unity 
(ideal target) at some definite angles and at some def- 
inite wavelengths. For example, TFCALC uses needle 
optimization to find the best thicknesses for such a merit 
function. We have preferred, however, to implement a 
gradient-based modified quasi-Newton algorithm (using 
the Fortran NAG libraries), for its consistency with the 
problem investigated, which is continuous. 

We have optimized the system [LH] N for N running 
from 5 to 20. When the thicknesses are kept equal 
by pairs, so as the structure retains its periodicity, we 
find thicknesses distributed around dh — 145 nm and 

„ — (N) 

dtf ~ 60 nm. The variation of 1Z with N is shown 

(TV) 

in Fig. 1, and can be numerically fitted to 1Z — 
0.774 + 0.365 x 10 -3 N; i.e., a linear dependence with an 
extremely small slope. We have also worked out the in- 
stance when all the thicknesses may vary independently 
(though now the system is not strictly periodic). The 
optimum thicknesses oscillate a lot, without a definite 
pattern (we do not list all of these values because of 

space limitation). However, as is clear from Fig. 1, 7Z 
shows an exponential increasing that can be suitably rep- 
resented by K {N) = 0.625 + 0.369[1 - exp(-0. 1327V)]. 

To gain further insights into these striking differences, 
in Fig. 2 we have plotted the reflectance VS 20 ^ for the 
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Fig. 2. (Color on line) Plot of the reflectance TZ^ 20 ^ as 
a function of the incidence angle 9 (in degrees) and the 
wavelength A (in fim). In a) we have the [LH] 20 system, 
while in b) the thicknesses are allowed to be different. 
At the bottom plane, we show the contour plots corre- 
sponding to a reflectance 0.99. The white lines delimit 
the zones in which this reflectance is greater than 0.99 
for all the angles of incidence. 



optimum thicknesses 2 ^ corresponding to N = 20 as a 
function of the angle of incidence and the wavelength A. 
At the bottom plane we have also included the contours 
of the regions where this reflectance is greater than 0.99. 
While for the [LH] 20 system, the top zone looks quite flat 
and very close to unity, the dips are very deep indeed. On 
the contrary, when we allow for different thicknesses, the 
top zone presents small ripples, but the dips are much 
less pronounced. Curiously enough, the [LH] 20 system 
gives a wider region lying above the 0.99 level. We have 
also marked two stripes (within parallel lines) in which 
the reflectance is greater than 0.99 for all the angles of 
incidence (and that, roughly speaking, could be identified 
with stop bands). Although of similar extension, they lie 
in quite different spectral ranges. We have repeated these 
calculations with other values of N, observing essentially 
the same kind of behavior. 

To proceed further, in Fig. 3 we have represented the 
angle-averaged reflectance [i.e., the magnitude in paren- 
theses in Eq. @] for the same two systems as in Fig. 2, in 
terms of the wavelength A. Obviously, the area under the 

curve is precisely 1Z^ 2 °\ It is clear that for different thick- 
nesses, this area is considerably bigger (so, it is really an 
optimum), while the [LH] 20 is better behaved in a nar- 
row range going from 0.57 to around 0.62 /im (which 
matches well with the stop band shown in Fig. 2. a). In 
other words, the optimum system remarkably improves 
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Fig. 3. (Color on line) Angle-averaged reflectance in 
terms of the wavelength A for the optimized [LH] 20 sys- 
tem (blue) and for the analogous system with different 
thicknesses (red). The line of reflectance equal to 0.95 is 
also shown. 

the behavior in the spectral wings, while it little affects 
(or even deteriorates) the behavior in the "good" central 
region. The ripples in both curves are caused by the dips 
appearing in Fig. 2 for each line of A constant. 

In this respect, we wish to note that, in our opinion, 
the notion of bandwidth becomes fuzzy for ODRs. Usu- 
ally^, it is defined as <5A = A+ — A_ [and sometimes 12 
normalized to the central wavelength (A++A_)/2], where 
A+ and A_ are the longer- and shorter-wavelength edges 
for given ODR bands [i.e., the two solutions of Eq. ([Tjl]. 
This is meaningful in the limit N — > oo, when these 
band edges make unambiguous sense, but fails for our 
more realistic situation. Other authors^ note that the 
common reflecting band for both polarizations and for 
angles up to a given 9m (which is defined by conven- 
tion) is [A_(#m), A+(0)] and the corresponding band- 
width <5A = A+(0) — A_(#m)- Again, this kind of defi- 
nition assumes the existence of a full band, which is only 
true for the strictly periodic case. 

We emphasize that all the relevant information about 
omnidirectionality is contained in the angle-averaged re- 
flectance. Furthermore, being an integral estimator, it 
does not rely on the values of the reflectance at some spe- 
cific relevant angles. For this reason, a sensible choice for 
defining the bandwidth is precisely the spectral range(s) 
for which this angle-averaged reflectance is bigger than a 
fixed threshold value. For example, if we set this value 
to, say 0.95, it is evident in Fig. 3 that the bandwidth 
for the [LH] 20 system is always poor. 

To sum up in a few words, we have exploited the no- 
tion of wavelength- and angle-averaged reflectance to ex- 
plore in a systematic way the performance of finite ODRs. 
Needless to say, our approach is general and can be ap- 



plied to other materials and other spectral regions. 
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